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FOREWORD 


This  interim  technical  report  was  prepared  by  Yale  University  Observatory, 
New  Haven,  Connecticut,  on  Contract  AF  33(657) -7271  for  the  Aeronautical  Re¬ 
search  Laboratories,  Office  of  Aerospace  Research,  United  States  Air  Force. 

The  research  reported  herein  was  accomplished  on  Task  7073-01,  "High  Energy 
Plasma  Generation  and  Control"  of  Project  7073,  "Research  on  Plasma  Dynamics" 
under  the  technical  cognizance  of  Dr.  Wolfgang  Braun  of  the  Plasma  Physics  Re¬ 
search  Laboratory  of  ARL. 


ii 


ABSTRACT 


This  report  describes  a  part  of  a  theoretical  investigation  of  the  theory  of 
the  emission  and  absorption  of  radiation  from  a  fully  ionized  plasma  in  the 
presence  of  a  magnetic  field.  Specifically,  the  bremsstrahlung  radiation  and 
departures  from  L.  T.  E.  are  discussed. 
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INTRODUCTION 


The  work  carried  out  under  this  contract  aims  at  developing  a  theory  of  emission 
and  absorption  of  radiation  in  fully  ionized  plasmas  the  presence  of  a  magnetic 
field.  The  pertinent  radiation  mechanisms  are  bremsstrahlung  and  cyclotron 
radiation.  A  general  theory  of  these  combined  mechanisms  must  take  into  account 
several  different  problems  which  determine  in  each  case  the  line  of  attack. 

At  present,  more  or  less  explicit  theories  are  available  for  special  cases  v'hich 
are  to  be  combined  and  expanded  to  cover  the  full  range  of  physical  possibilities. 

It  was  decided  to  attack  the  following  problems  in  turn: 

1.  Emission  and  absorption  of  brems strahlung  following  Coulomb  interactions  is 
influenced  by  the  presence  of  a  magnetic  field,  whereas  at  the  same  time  the 
cyclotron  line  is  broadened  by  the  Coulomb  interactions.  A  theory  of  the  combined 
effects  was  to  be  developed,  neglecting  for  the  time  being  refractive  indices  differ¬ 
ing  from  unity,  departures  from  thermodynamic  equilibrium  and  relativistic 
corrections.  2.  The  influence  of  the  refractive  index  was  to  be  considered.  This 
effect  is  a  manifestation  of  the  collective  behavior  of  the  plasma.  Since  a  systematic 
theory  is  still  missing,  it  was  judged  most  advisable  to  begin  work  with  simplest 
case,  namely,  the  influence  of  shielding  corrections  on  the  bremsstrahlung  cross 
section.  3.  Departures  from  thermodynamic  equilibrium  do  not  affect  radiative 
interactions  of  single  particles,  but  only  emission  and  absorption  of  the  plasma 
as  a  whole.  Since  this  is  essentially  a  transfer  problem,  it  could  be  separated 
from  the  remainder  of  the  analysis.  4.  Relativistic  corrections  in  a  general 
manner  are  very  cumbersome  to  treat.  It  was  therefore,  decided  to  treat 
relativistic  corrections  in  each  of  the  general  cases  under  study  separately,  and 
at  a  later  stage  of  the  investigation. 
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In  the  period  covered  by  this  report,  a  general  theory  of  brems strahlung 
and  cyclotron  radiation  for  unity  refractive  indices  was  developed  and  is  briefly 
described  in  Section  1.  The  theory  of  brems  strahlung  emission  was  expanded 
to  include  refractive  indices  differing  from  unity,  leading  to  extensive  numerical 
work  which  is  now  under  way;  c.  f.  Section  2.  Departures  from  thermodynamic 
equilibrium  were  treated  in  a  very  general  manner,  and  an  exhaustive  formula¬ 
tion  for  continuous  energy  developed.  The  results  are  summarized  in  Section  3. 
1.  RADIATIVE  COULOMB  INTERACTION  IN  THE  PRESENCE  OF  A  MAGNETIC 
FIELD. 

The  purpose  of  this  study  is  to  obtain  a  theory  of  Coulomb  interactions  in  the 
presence  of  a  magnetic  field.  The  essential  parameter  is  the  collision  cross 
section  of  the  free  electrons  with  ions  in  the  magnetic  field.  The  method  to  be 
followed  is  to  compute  from  the  wave  functions  of  free  electrons  in  a  magnetic 
field  the  matrix  elements  for  dipole  radiation,  where  the  Coulomb  interactions 
are  treated  as  perturbations  of  the  time  varying  part  and  of  the  space  varying 
part  of  the  wave  functions  separately. 

The  statistical  superposition  of  these  perturbations  is  carried  out  in  the 
manner  given  by  Lindholm^  and  leads  directly  to  the  expressions  for  line  shift 
and  half  width.  The  two  effects  are  separable.  If  one  neglects  the  effect  of  the 
Coulomb  interactions  on  the  space  depending  part  of  the  wave  functions  (adia¬ 
batic  case  in  the  impact  theory  of  line  broadening),  only  a  certain  fraction  of  the 
cyclotron  line's  half  width  and,  for  that  matter,  of  the  scattering  cross  section 
of  free  electrons  in  a  magnetic  field,  is  obtained.  However,  a  line  shift  is  pre¬ 
dicted  which  does  not  follow  from  the  non-adiabatic  case,  i.  e.  ,  the  effect  of  the 
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Coulomb  interactions  on  the  space  depending  part  of  the  wave  functions. 


In  the  following,  the  adiabatic  treatment  is  reviewed  briefly  for  which 
definitive  results  are  now  available.  We  mention,  however,  that  preliminary 
results  indicate  that  the  total  cross  section  following  from  the  non-adiabatic 
treatment  does  not  depend  critically  on  the  magnetic  field,  since  in  most  cases 
the  dependence  on  the  magnetic  field  strength  is  logarithmic.  Detailed 
computations  covering  the  general  case  are  scheduled  for  the  following  research 
period. 
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The  spatial  part  of  the  wave  functions  reads 

(n!s!r''*e-'’'><"-s)/2  Q”-^(p). 

The  normalization  length  L  is  arbitrary  and  does  not  appear  in  the  final  result. 
The  parameter 

X=eH/2c1i 


contains  the  magnetic  field.  yO  is  defined  by 

p-.yx^ 


(3) 


r  being  the  radius  vector,  k  is  a  continuous  variable,  s  is  integer  and  varies 
from  0  to  ,  n-s  varies  from  -«>  to  n. 
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is  the  generalized  Laguerre  polynomial. 

So  far,  we  have  considered  the  wave  function  of  electrons  subject  only  to  a 
constant  magnetic  field.  We  now  proceed  to  include  the  effects  of  interactions 
between  the  electron  and  positive  ions  in  the  two  particles  approximation.  The 
interaction  adds  a  time  and  space  varying  part  H(r,t)  to  the  Hamiltonian  in  the 
Schrodinger  equation  which  now  reads 


+  H(Ti) 


^  =  -iti 


at 


(5) 


In  Eq  (5),  ijj  is  the  wave  function  of  the  electron,  p  the  momentum  operator.  The 
vector  potential  A  and  the  magnetic  field  H  are  related  by 


Ax  A  =7?  (6) 

Assuming  that  at  the  initial  time  there  are  no  perturbations  present,  we  have  as 
solution  the  value  of  4*  from  Eq  (1).  The  time  varying  part  can  be  written  in  the 
form  exp  [iEt/h]. 

The  problem  of  solving  Eq  (5)  is  exceedingly  complex.  The  following  simpli- 
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fications  were,  therefore,  introduced  ,  Firstly,  we  represent  the  entire  set  of 
possible  electron-ion  interactions  by  the  same  electron  wave  function,  but  by 
varying  ion  positions.  The  wave  function  of  the  electron  can  then  be  chosen  as 
a  function  which  is  sharply  peaked  at  the  Larmor  radius.  Mathematically,  this 
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implies  that 


S  =  0  (7) 

or,  in  physical  terms,  that  the  origin  is  chosen  coincident  with  the  center  of  the 
orbit. 

Secondly,  as  already  mentioned  above,  we  assume  that  the  effect  of  the  elec¬ 
tron  ion  interaction  is  to  change  the  time  varying  part  of  the  wave  functions  only 
(adiabatic  case).  The  classical  counterpart  is  to  neglect  a  spatial  change  in  the 
orbit  of  the  particle  due  to  the  interaction.  If  there  were  no  magnetic  fields 
present,  this  orbit  would  be  a  straight  line. 

Lindholm^  has  developed  a  theory  for  the  broadening  of  spectral  lines  in  the 
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case  of  interactions  producing  a  "random  phase  change"  .  This  type  of  interac¬ 
tion  produces  a  certain  change  in  the  eigenfrequency  of  the  oscillator,  whereby 
in  a  large  number  of  interactions  the  distribution  of  frequency  modulations  is 
random.  In  essence,  the  postulated  randomness  amounts  to  saying  that  the  inter¬ 
actions  are  uncorrelated,  and  that  their  effects,  therefore,  are  additive.  It  is 
not  the  scope  of  this  report  to  review  the  several  steps  leading  to  Lindholm's 
final  result  which  reads  for  the  spectrum  of  the  cyclotron  line 


r  21“' 

®  K  +((*<-0)8 -Uj)  J  doj 


(8) 
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The  line  contour  represented  by  Eq  (8)  is  very  close  to  a  Lorentzian  form,  except 
that  in  addition  to  the  classical  broadening  represented  by  the  parameter 
there  occurs  a  shift  of  the  resonance  frequency 


cu^  =  eH/mc 
c 


(9) 


which  is  given  by  u^. 

Half-width  Uj^  and  shift  u^  were  calculated  with  the  aid  of  the  wave  functions 
of  Eq  (1)  and  the  Coulomb  perturbations  in  Schrodinger '  s  Eq  (5),  subject  to  the 
simplifications  outlined  above.  The  results  are  given  in  Table  1  and 

Figure  1  (u^).  Whereas,  the  total  cross  section  will  differ  from  u^  by  the  non- 

adiabatic  contribution,  the  line  shift  is  the  definitive  result.  For  a  combination 

10-3  3 

of  i oft  number  density  n.  =  10  cm  and  H=  10  Gauss,  the  line  shift  is  about 

1 

3  percent  of  the  gyrofrequency.  Hence,  it  should  be  observable  under  suitable 
experimental  conditions. 

2.  BREMSSTRAHLUNG  CROSS  SECTIONS  FOR  FREQUENCIES  CLOSE  TO  THE 
'  PLASMA  FREQUENCY 

The  conventional  theory  of  bremsstrahlung  interactions  assumes  that  the 
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refractive  index  (ri)  is  effectively  unity  .  Since  (in  the  absence  of  a  magnetic 
field),  the  r.  i.  (I,  Eq  77) 

n=  [l-ajp/cu^]'^^  (10) 

I 

i 

I 

i 

j  where  to^  is  the  plasma  frequency,  the  condition  of  r.  i.  unity  requires 
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(11) 


(0  » (Jjp 

On  the  other  hand,  the  inequality  (11)  is  required  in  order  to  justify  the  use 
of  a  pure  Coulomb  potential  for  the  electron  ion  interaction  without  shielding 
corrections.  In  fact,  the  departures  of  the  r.  i.  from  unity  are  due  to  the  shield¬ 
ing  effects  of  the  plasma  electrons  on  the  ion's  potential^.  Hence,  accounting 
for  values  of  n^l  means  to  include  shielding  effects  in  the  calculation  of 
bremsstrahlung  cross  sections. 

In  first  approximation,  the  shielding  effects  can  be  included  by  cutting  off 
the  Coulomb  potential  at  the  Debye  distance^,  and  by  using  this  cut-off  potential 
instead  of  the  pure  Coulomb  potential  in  the  conventional  classical  calculation  of 
the  bremsstrahlung  cross  sections. 

This  calculation  must  then  be  altered  in  two  ways.  Firstly,  the  time  integra¬ 
tion  in  the  evaluation  of  the  Fourier  components  I,  Eqs  (8)  and  (9)  must  be  terminated 
at  the  finite  times  when  the  electron  enters  and  leaves  the  Debye  sphere, 
respectively.  Secondly,  only  interactions  of  electrons  passing  the  ion  at  distances 
(collision  parameters)  less  than  the  Debye  distance  are  counted  for  the  average 
cross  section. 

Making  use  of  the  formulation  developed  in  I,  we  have  in  the  straight  line 
approximation^.  Section  4]  which  is  valid  throughout  the  r.f.  region  and,  thus, 
in  the  neighborhood  of  the  plasma  frequency, 

2  "^^0 

x(a))=  ■  ^  h  /  cos(^leu)  [l  +  u^r’  du  (12) 

°  -Uo 
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and 


where 


TTmvQb 


/sin(il€u)[l+u^]  udu, 

-Uo 


1/2 


(14) 


e,  the  eccentricity  of  the  hyperbolic  orbit,  and  the  collision  parameter  b  are 
related  for  a  given  electron  velocity  v^  by 


rnnvo*  i 

2 

‘  bcu 

L  Ze"  J 

(15) 


follows  from  €  by  inserting  the  Debye  distance 


b„=(kT)''V2Tr''*eN. 


(16) 


[c.f.  I,  Eq  70]  for  b.  Since  ^  ^  1, 


b 


(17) 


The  integrals  (12  and  13)  were  computed  numerically.  The  spectrum  coefficient 


(18) 
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[I,  Eq  71]  must  then  be  averaged  with  respect  to  all  collision  parameters  or, 
what  amounts  to  the  same,  with  respect  to  the  eccentricities  e.  The  pertinent 
average  reads  [I,  Eq  59] 

m 

TT  <Q^)hdcu  =  27r  Nj  VqTT  bdb'daj  (i9) 

0 

3 

Here,  N.  is  the  number  of  ions  per  cm  .  Note  the  finite  upper  limit  due  to  the 
shielding  of  the  ion  potential. 

Preliminary  results  for  a  few  test  cases  show  that  the  average  spectrum 

coefficient  is  a  smooth  curve  decreasing  rapidly  for  frequencies  above  the  plasma 

frequency.  This  behavior  is  expected  on  general  grounds  and  agrees,  at  least 
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qualitatively,  with  the  results  of  Dawson  and  Oberman  . 

3.  DEPARTURES  FROM  LTE 

The  importance  of  departures  from  local  thermodynamic  equilibrium  (LTE) 

was  first  recognized  in  connection  with  certain  astr ophysical  problems,  in 

particular,  in  connection  with  the  interpretation  of  the  spectrum  of  the  outer 
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solar  atmosphere  ,  Stellar  plasmas  are,  in  the  observable  layers,  rather  weakly 
ionized,  so  that  the  major  part  of  the  radiation  spectrum  consists  of  spectral 
lines,  i.  e.  ,  emissions  and  absorptions  involving  only  bound  states.  Departures 
from  LTE  thenaffect  primarily  these  bound  states,  with  the  free  energy  states  in 
LTE.  In  physical  terms,  this  statement  implies  that  the  free  energy  states  are 
distributed  according  to  a  Maxwell  distribution,  whereas,  emission  and  absorp¬ 


tion  are  connected  by  Kirchhoff’s  law 


^U)  ~ 


(20) 


e  is  the  emission  per  cm  ,  sec.  and  steradian,  k  is  the  absorption  coefficient 
O)  u> 

per  cm,  B  is  Planck’s  fxinction. 

CO 

In  many  domaines  of  laboratory  physics,  however,  plasmas  are  completely 
ionized,  and  the  population  of  bound  states  can  be  altogether  neglected.  Then, 
only  the  free  electrons  are  involved  in  radiative  interactions,  and  departures 
from  LTE  in  the  radiation  fi^ld  have  their  counterpart  in  departures  from  a 
Maxwellian  distribution  of  the  free  energy  states. 

The  purpose  of  this  outline  is  to  give  a  short  description  of  a  general  formu¬ 
lation  of  departures  from  LTE  in  free-free  emission  and  absorption.  We  restrict 
our  considerations  to  steady  states,  where  the  change  in  the  population  of  states 
is  (during  observational  times)  small  compared  with  the  populations  themselves. 
The  effects  of  non-radiative  interactions  are  of  minor  importance  for  the 
problems  treated  presently,  since  the  conclusions  of  this  Section  do  not  depend 
on  the  presence  or  absence  of  non-radiative  interactions.  A  steady  state,  on  the 
other  hand,  is  the  case  of  greatest  interest  under  laboratory  conditions. 

Radiative  interactions  affect  the  particle  distribution  and  the  spectrum  of  the 
radiation  field.  There  are,  therefore,  always  two  sets  of  conditions  to  be 
considered,  derived  from  the  behavior  of  the  radiation  field,  and  from  the  be¬ 
havior  of  the  particles.  Neglecting  one  set  of  conditions  does  mathematically  not 

lead  to  any  contradiction,  but  leads  to  physically  invalid  solutions.  This  fact  has 
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been  pointed  out  previously  . 

The  steady  state  condition  for  the  particle  states  implies  that  the  number  of 
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transitions  leading  into  a  given  state  of  energy  from  all  states  above  (Ej.) 

and  below  (E  )  must  equal  the  number  of  transitions  out  of  E  .  Using  the 

o 
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analysis  published  previously  ,  we  find  the  steady  state  condition 


/dE  QiE  + 

•'En=  const 

fo 

/dE,Q{Eo,ii){Eo,I2.} 


-0  Et 


(21) 


=  0 


where 


Eo=  const 


=  Ej-Eo  =Eo“^2 


(22) 


The  bracket  symbol  is  defined  by 


{E.fl}  =  N(E)  +  •  4^^[n(E)-  N(E-il)]  ,33, 

which  contains  spontaneous  and  induced  emissions,  and  absorptions,  already 
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reduced  to  the  probability  of  spontaneous  emission 

Q(E,il)  dEd^i  [sec"'  sterad"'] .  ^24) 

2 

I  d^2  is  the  specific  intensity  of  the  impinging  radiation  field  (erg /cm  sec 

CO 

steradian),  N(E)  is  the  (arbitrary)  number  distribution  of  electrons  per  energy 
interval  and  cm^. 

If  no  additional  conditions  are  imposed  on  the  system  (see  below),  the 
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consideration  of  the  radiation  fields  adds  the  conservation  law  of  energy  as  a 
kind  of  boundary  condition.  In  the  above  formalism,  this  condition  reads 


0  = 


00  00 

/dEo  /dE 

°  Eo 


Q(E|  ,fi){E|,^  =  /dEo  /dEj  •  Q(Eo,il){Eo,^i} 


(25) 


The  set  of  Eqs  (21)  and  (25)  has  the  trivial  solution 


{E.a}sO 


for  all 


E, 


(26) 


It  can  be  easily  shown  that  Eq  (26)  is  equivalent  to  LTE,  i.  e.  ,  to  a  system  where 

N(E)  is  Maxwellian,  I  is  Planckian, 

w 

It  can  also  be  shown  with  a  little  algebra  that  in  a  system  which  is  optically 
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thick  in  the  whole  spectrum,  i.  e.  ,  in  a  system  where 


(27) 


(S  is  the  source  function),  everywhere,  the  conservation  law  of  total  radiative 

CO 

energy  Eq  (25)  is  replaced  by  the  much  stronger  condition  that 


The  only  general  solution  of  Eq  (2  8)  is  the  trivial  solution  (26)  for  LTE. 
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Hence,  we  conclude  that  a  plasma  which  is  optically  thick  in  the  whole  continuum 
cannot  depart  from  LTE  under  stationary  conditions.  Since  Eq  (28)  does  not 
depend  on  the  presence  of  non-radiative  transitions,  the  above  stated  conclusion 
is  valid  generally. 

At  present,  systems  are  investigated  in  which  Eq  (27)  holds  only  in  certain 
regions  of  the  spectrum. 
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